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I. INTRODUCTION 

The determination of the properties of nuclear matter 
as function of density and temperature is a fundamen- 
tal task in nuclear and subnuclear physics. Heavy ion 
collisions experiments open the possibility to investigate 
strongly interacting compressed nuclear matter exploring 
in laboratory the structure of the QCD phase diagram 0- 
4] . The extraction of information about the Equation of 
State (EOS) at different densities and temperatures by 
means of intermediate and high energy heavy ion colli- 
sions is a very difficult task and can be realized only indi- 
rectly by comparing the experimental data with different 
theoretical models, such as, for example, fluid-dynamical 
models. The EOS at density below the saturation den- 
sity of nuclear matter (po m 0.16 fm~ 3 ) is relatively well 
known due to the large amount of experimental nuclear 
data available. At larger density there are many uncer- 
tainties; the strong repulsion at short distances of nuclear 
force makes, in fact, the compression of nuclear matter 
quite difficult. However, in relativistic heavy ion colli- 
sions the baryon density can reach values of a few times 
the saturation nuclear density and/or high temperatures. 
The future CBM (Compressed Baryonic Matter) experi- 
ment of FAIR (Facility of Antiproton and Ion Research) 
project at GSI Darmstadt, will make it possible to cre- 
ate compressed baryonic matter with a high net baryon 
density [5|-|7[. In this direction interesting results have 
been obtained at low SPS energy and are foreseen at a 
low-energy scan at RHIC [8rfl2l|. 

On the other hand, the information coming from ex- 
periments with heavy ions at intermediate and high en- 
ergy collisions is that the EOS depends on the energy 
beam but also on the electric charge fraction Z/A of the 
colliding nuclei, especially at not too high temperature 
[HI, [lij]. Moreover, the analysis of observations of neu- 
tron stars, which are composed of /3-stable matter for 
which Z/A < 0.1, can also provide hints on the structure 
of extremely asymmetric matter at high density [ill (3 ■ 

To well understand the structure of the phase diagram 
and the supposed deconfinement quark-gluon phase tran- 
sition at large density and finite temperature, it is crucial 
to know accurately the EOS of the hadronic as well as 



the quark-gluon phase. Concerning the hadronic phase, 
hadron resonance gas models turned out to be very suc- 
cessful in describing particle abundances produced in (ul- 
tra)relativistic heavy ion collisions |17Hl9| . In this frame- 
work, to take phenomenologically into account the in- 
teraction between hadrons at finite densities, finite size 
corrections have been considered in the excluded volume 
approximation [2p|-|25j|. 

From a more microscopic point of view, the hadronic 
EOS should reproduce properties of equilibrium nuclear 
matter such as, for example, saturation density, bind- 
ing energy, symmetric energy coefficient, compression 
modulus. Some other constraints on the behavior of 
the EOS come out from analysis of the experimental 
flow data of heavy ion collisions at intermediate en- 
ergy [2t| and, moreover, there are different indirect con- 
traints/indications related to astrophysical bounds on 
high density /3-equilibrium compact stars [HI, [27|. In 
connection with these matters, Walecka type Relativistic 
Mean-Field (RMF) models have been widely successfully 
used for describing the properties of finite nuclei as well 
as dense and finite temperature nuclear matter [28rl32| . 
It is relevant to point out that such RMF models usually 
do not respect chiral symmetry. Furthermore, the repul- 
sive vector field is proportional to the net baryon density, 
therefore, standard RMF models do not appear, in prin- 
ciple, fully appropriate for very low density and high tem- 
perature regime. In this context, let us observe that has 
been recently proposed a phenomenological RMF model 
in order to calculate the EOS of hadronic matter in a 
broad density-temperature region by considering masses 
and coupling constants depending on the cr-meson field 
[33| . In that approach, motivated by the Brown- Rho scal- 
ing hypothesis, a not chiral symmetric model simulates a 
chiral symmetric restoration with a temperature increase. 
On the other hand, sophisticated relativistic chiral SU(3) 
models are also developed to take into account particle 
ratios at RHIC and baryon resonances impact on the chi- 
ral phase transition [34], [35J . 

In regime of finite values of density and temperature, 
a state of high density resonance matter may be formed 
and the A(1232)-isobars degrees of freedom are expected 
to play a central role [3(|. Transport model calculations 
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and experimental results indicate that an excited state 
of baryonic matter is dominated by the A-resonance at 
the energy from AGS to RHIC [I3-E1- Moreover, in the 
framework of the non-linear Walecka model, it has been 
predicted that a phase transition from nucleonic mat- 
ter to A-exited nuclear matter can take place and the 
occurrence of this transition sensibly depends on the A- 
mesons coupling constants [43|, |44| . Referring to QCD 
finite-density sum rules results, which predict that there 
is a larger net attraction for a A-isobar than for a nucleon 
in the nuclear medium [45j , the range of values for the 
A-mesons coupling constants has been confined within a 
triangle relation [46|. Whether stable A-exited nuclear 
matter exists or not is still a controversial issue because 
little is actually known about the A coupling constants 
with the scalar and vector mesons. In any case, it has 
been pointed out that the existence of degrees of free- 
dom related to A-isobars can be very relevant in rela- 
tivistic heavy ion collisions and in the core of neutron 
stars 0, H3 |48|. Although several papers have inves- 
tigated the influence of A-isobars in the nuclear EOS, 
we believe that, especially in presence of asymmetric and 
strange hadronic matter, a systematic investigation at 
finite densities and temperatures has been lacking. 

In this paper we are going to study the hadronic EOS 
by means of an effective RMF model with the inclusion of 
the full octet of baryons, the A-isobars degrees of freedom 
and the lightest pseudoscalar and vector mesons. These 
last particles are considered in the so-called one-body 
contribution, taking into account their effective chemi- 
cal potentials depending on the self-consistent interaction 
between baryons. The main goal is to investigate how the 
constraints on the global conservation of the baryon num- 
ber, electric charge fraction and strangeness neutrality, in 
the presence of A-isobars degrees of freedom, hyperons 
and strange mesons, influence the behavior of the EOS in 
regime of finite values of baryon density and temperature. 
Moreover, we plan to show the relevance of A-isobars for 
different coupling constants and how its presence influ- 
ences several particle ratios and strangeness production 
for three different parameters sets, compatible with ex- 
perimental constraints. 

The paper is organized as follows. In Section II, 
we present the model with a detailed discussion on 
the hyperons-mesons couplings and on the chemical 
equilibrium conditions. In order to better clarify the 
role of A-isobars and strange particles in symmetric and 
asymmetric nuclear matter, our results are presented 
in Section III which is divided in three subsections: 
in A) we study the equation of state of nucleons and 
A-isobars in symmetric nuclear matter at zero and 



finite temperature; in B) we extend the investigation 
by including hyperons, non-strange and strange mesons 
in asymmetric nuclear matter and by requiring the 
zero net strangeness condition. Strangeness production 
and different particle-antiparticle ratios are considered 
in subsection C). Finally, we summarize our main 
conclusions in Section IV. 



II. HADRONIC EQUATION OF STATE 

The basic idea of the RMF model, first introduce d by 
Walecka and Boguta-Bodmer in the mid-1970s [HI, l50f . 
is the interaction between baryons through the exchange 
of mesons. In the original version we have an isoscalar- 
scalar a meson field which produces the medium range 
attraction and the exchange of isoscalar- vector ui mesons 
responsible for the short range repulsion. The satura- 
tion density and binding energy per nucleon of nuclear 
matter can be fitted exactly in the simplest version of 
this model but other properties of nuclear matter, as e.g. 
incompressibility, cannot be well reproduced. To over- 
come these difficulties, the model has been modified in- 
troducing in the Lagrangian two terms of self-interaction 
for the a mesons which are crucial to reproduce the em- 
pirical incompressibility of nuclear matter and the effec- 
tive mass of nucleons M^. Moreover, the introduction 
of an isovector-vector p meson allows to reproduce the 
correct value of the empirical symmetry energy [5l| and 
an isovector-scalar field, a virtual ao(980) (5-meson, has 
been been studied for asymmetric nuclear matter and for 
heavy ion collisions [l3|, [52[ . 

The total Lagrangian density C can be written as 

£ = £om + -Ca + £qfm , (1) 

where £ om stands for the full octet of lightest baryons 
(p, n, A, S + , S°, £ _ , 5°, S _ ) interacting with a, ui, p, S 
mesons fields; £a corresponds to the degrees of freedom 
for the A-isobars (A ++ , A + , A , A~) and £ q f m is related 
to a (quasi)free gas of the lightest pseudoscalar and vector 
mesons with an effective chemical potential (see below 
for details). In the regime of density and temperature we 
are mostly interested, we expect that the inclusion of the 
other decuplet baryons will produce only a small change 
in the overall results. 

The RMF model for the self-interacting full octet of 
baryons (J p — l/2 + ) was originally studied by Glenden- 
ning [EH with the following standard Lagrangian 
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£om = ^2 [i 7p ^ - ( M k ~ gakcr ~ gsk t- 6) - g uk 7 p w M - g pk j p f • p M ] ^> fe + - (d^crd^a - to 2 er 2 )-U{a) 
k 

+~ ml + \ c (gl N c^) 2 + \ m 2 p M • + ~(djd»6- mj P) - \g^& v , (2) 

I 



where the sum runs over all baryons octet, M k is the 
vacuum baryon mass of index k, the quantity t denotes 
the isospin operator which acts on the baryon and the 
field strength tensors for the vector mesons are given by 
the usual expressions F^ v = d p io v — d v uj p , G pv = d^p v — 
dvpu- The U{a) is a nonlinear self-interaction potential 
of a meson 



u (c) = (g aN cr) 3 + i b (g aN a 4 ) , 



(3) 



lowing form 

(i-t^ - M* k - g^u - g pk j°t 3k p)^ k = , (5) 

(i^d» - Ml - 5w A7°w) = , (6) 

ml<r + agl N a 2 +bg* N cr 3 =^2g at pf , (7) 

i 

mluj + c g* N uj 3 = g^ pf , (8) 

i 

= J2g P it 3i pf, (9) 

(10) 



m 2 p- 



mgS = Y gsi hi pf , 



introduced by Boguta and Bodmer [50( to achieve a rea- 
sonable compressibility for equilibrium normal nuclear 
matter. We have also taken into account the additional 
self-interaction w meson field: c (g 2 N w Al a; M ) 2 /4 suggested 
by Bodmer (54j to get a good agreement with Dirac- 
Briickner calculations at high density and to achieve a 
more satisfactory description of the properties of finite 
nuclei in the mean field approximation. 

By taking into account only the on-shell As, the La- 
grangian density concerning the A-isobars can be ex- 
pressed as [il] 



Ca=iP Av Vl^ ~ ( m a - g*Acr) - gu,Al^W A ) ( 4 ) 



where ipl is the Rarita-Schwinger spinor for the A 
baryon. Due to the uncertainty on the A-meson cou- 
pling constants, we limit ourselves to consider only the 
coupling with the a and u> mesons fields, more explored 
in the literature (13-111. 



In the RMF approach, baryons are considered as Dirac 
quasiparticles moving in classical mesons fields and the 
field operators are replaced by their expectation values. 
In this context, it is relevant to remember that the RMF 
model does not respect chiral symmetry and the contri- 
bution coming from the Dirac-sea and quantum fluctua- 
tion of the meson fields are neglected. As a consequence, 
the field equations in RMF approximation have the fol- 



where a = (a), uj = (uj°), p — (p®) and 5 = {S3) are the 
nonvanishing expectation values of mesons fields. The 
effective mass of fc-th baryon octet, comparing in Eq.([5]), 
is given by 



Ml=M k - g ak a - gs k t 3k 5 , 



(11) 



and the effective mass of A-isobar, comparing in Eq. (JS|) , 
is given by 



Ml = M A - gaA<J • 



(12) 



In the mesons fields equations, Eqs.([T))- (fTU)) , the sums run 
over all considered baryons (octet and As) and pf and 
pf are the baryon density and the baryon scalar density 
of the particle of index i, respectively. They are given by 



Pi = li 



Pf = 7* 



d 3 fc 



(2tt) 3 
d 3 fc M, 



[fi(k)-fi(k)}, 



(13) 



(2tt)3 E* 



[fm+fm, (14) 



where ji = 2 Ji + l is the degeneracy spin factor of the i-th 
baryon (7 oc tct = 2 for the baryon octet and 7a = 4) and 
fi(k) and fi(k) are the fermion particle and antiparticlc 
distributions 



fi(k) : 

7i(k) 



1 



exp{(E*(k)-p*)/T} + l 
1 



(15) 
(16) 



exp{(£?(fc) + M?)/T} + l 

The baryon effective energy is defined as E*(k) = 
k 2 + Mi* 2 . The chemical potentials p,i are given in 
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terms of the effective chemical potentials p* by means of 
the following relation 



Mi = Pi + 9ui u + g pi hi p , 



(17) 



where is the third component of the isospin of i-th 
baryon. 

Because we are going to describe the nuclear EOS at 
finite density and temperature with respect to strong in- 
teraction, we have to require the conservation of three 
"charges": baryon number (B) , electric charge (C) and 
strangeness number (S). Each conserved charge has a 
conjugated chemical potential and the systems is de- 
scribed by three independent chemical potentials: pb, 
pc and Ms- Therefore, the chemical potential of particle 
of index i can be written as 



Hi = hi p B + c t p c + si ps ■ 



(18) 



where bi, Ci and Sj are, respectively, the baryon, the elec- 
tric charge and the strangeness quantum numbers of i-th 
hadronic species. 

The thcrmodynamical quantities can be obtained from 
the grand potential tts in the standard way. More explic- 
itly, the baryon pressure Pb = — ^b/V and the energy 
density cb can be written as 



Pr = 



IE 



3 

% 

1 2 2 



d 3 /c fc 2 



(2tt)3 Ef(k) 



[fi(k) + fi(k)} 



U{<J) + i ml uj' 2 + i c (g uN u f 



1 2 2 1 2 c-2 
m p P ~ o 6 ' 



(19) 



£B = 



£7. 

i 

1 



1 2 2,1 2 2 

2 m P P + 2 m s S 



(20) 



The numerical evaluation of the above thermodynam- 
ical quantities can be performed if the meson-nucleon, 
meson-A and meson-hyperon coupling constants are 
known. Concerning the meson-nucleon coupling con- 
stants (gaN, g^N, gpN, gsN), they are determined to re- 
produce properties of equilibrium nuclear matter such as 
the saturation densities, the binding energy, the symmet- 
ric energy coefficient, the compression modulus and the 
effective Dirac mass at saturation. Due to a valuable 
range of uncertainty in the empirical values that must be 
fitted, especially for the compression modulus and for the 
effective Dirac mass, in literature there are different sets 
of coupling constants. In Table HI we report the parame- 
ters sets used in this work. The set marked GM3 is from 
Glendenning and Moszkowski (55[, that labelled NLpJ is 
from Ref.fllii] and TM1 from Ref.[57|. As we will see 
in next Section, we have limited our investigation to these 
three parameters sets that are compatible with interme- 
diate heavy ion collisions constraints and extensively used 



in various high density astrophysical applications. Let us 
remark that the first two parameters sets have the same 
saturated compressibility K and an almost equal value of 
the nucleon effective mass M^, significantly larger than 
the TM1 one. Therefore, the GM3 and NLpS models will 
fail to reproduce the correct spin-orbit splittings in finite 
nuclei [58|. On the other hand, the TM1 parameters set 
have a larger value of K but a sensibly lower value of 
M^. As we will see, these different saturation properties 
of nuclear matter are strongly correlated to the formation 
of A-isobar matter at finite density and temperature. 

The implementation of hyperons degrees of freedom 
comes from determination of the corresponding meson- 
hyperon coupling constants that have been fitted to hy- 
pernuclear properties. Following Ref.s[59-63j, we can use 
the SU(6) simple quark model and obtain the relations 



1 



1 



1 



9pN = 2-9pS = 5pH , g P A = , 
1 

9sn = 2#<5s = gs~. , gsA = . 



(21) 



In addition, we can fix the scalar a meson-hyperon {g a y) 
coupling constants to the potential depth of the corre- 
sponding hyperon in normal dense matter taking into 
account the following recent results [63l466| 



C/^ = -28MeV, 



u£=- 



30MeV, U^ = -18MeV. (22) 



In Table HH we report the obtained ratios x a y — 
g<rY I g<rN and the vacuum hyperon masses are listened 
in Table IIII1 In this context, let us mention that the 
two additional mesons fields /o(975) and </>(1020), usu- 
ally introduced to simulate the hyperon- hyperon attrac- 
tion observed in A — A hypernuclei (62l . [63| , do not play 
a significant role in the considered range of density and 
temperature and, therefore, we have neglected their con- 
tributions. 

As discussed in the Introduction, the aim of this work 
is to describe the EOS at finite values of density and tem- 
perature. Especially at low baryon density and high tem- 
perature, the contribution of the lightest pseudoscalar 
and vector mesons to the total thermodynamical po- 
tential (and, consequently, to the other thermodynam- 
ical quantities) becomes very important. On the other 
hand, the contribution of the 7r mesons (and other pseu- 
doscalar and pseudovector fields) vanishes at the mean 
field level. From a phenomenological point of view, we 
can take into account the meson particles degrees of free- 
dom by adding their one-body contribution to the ther- 
modynamical potential, i.e. the contribution of an ideal 
Bose gas with an effective chemical potential /x* , depend- 
ing self-consistently from the mesons fields. Following 
this scheme, we can evaluate the pressure Pm, the en- 
ergy density £m and the particle density pf 1 of mesons 
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TABLE I: Nuclear matter properties and nucleon coupling constants of the parameters sets used in the calculation. The energy 
per particle is _E/A=16.3 MeV, calculated at the saturation density po with a compression modulus K and effective mass 
(the nucleon mass Afjv is fixed to 939 MeV for GM3 and NLp<5, Afjv = 938 MeV in TM1 parameters set). The symmetry energy 
is denoted by a aym . In the parameters set NL/9<5 the additional coupling to S meson is fixed to gsN/ms =3.162 fm. 





Po 


K 


M* N /M N 


Q>sym 


9(7 N 
ma 


9u,N 


9pN 

m p 


ffl 


b 


c 




(fm" 3 ) 


(MeV) 




(MeV) 


(fm) 


(fm) 


(fm) 


(fm- 1 ) 






GM3 


0.153 


240 


0.78 


32.5 


3.151 


2.195 


2.189 


0.04121 


-0.00242 




NLp<5 


0.160 


240 


0.75 


30.5 


3.214 


2.328 


3.550 


0.0330 


-0.0048 




TM1 


0.145 


281 


0.63 


36.9 


3.871 


3.178 


2.374 


0.00717 


0.00006 


0.00282 



TABLE II: Ratios of the scalar er-meson coupling constants 
for hyperons: x„y = g^y/g^N- 











GM3 


0.606 


0.328 


0.322 


NLp<5 


0.606 


0.361 


0.320 


TM1 


0.616 


0.447 


0.319 



tive chemical potential can be written as 

Kr+(p+) = V*p ~ 

= fi p - p n - g pN p , 



(27) 



where the last equivalence follows from Eq. (fT7|) . 
Analogously, by setting i wA = 9u>a/ g^N, we have 



P 



M 



It. 



d 3 fc 



^K+(K' + ) =/ X p — MA(EO) 



gj(k). 



(2tt) 3 Ej(k) 
j 

EM = 22l] J 7^A3 E ii k ) 9A k ) > 

d 3 /c 



Pf = 7, 



(2tt)* 



g 3 {k). 



(23) 
(24) 
(25) 



where jj = 2 Jj + 1 is the degeneracy spin factor of the 
j-th meson (7=1 for pseudoscalar mesons and 7 = 3 
for vector mesons), the sum runs over the lightest pseu- 
doscalar mesons (tt, K, K, 77, rj') and the lightest vector 
mesons (p, oj, K*, K , <f>), considering the contribution 
of particle and antiparticle separately. In Eq.s (|2"31 - ([23)1 
the function gj(k) is the boson particle distribution (the 
corresponding antiparticle distribution g~j(k) will be ob- 
tained with the substitution p* — ¥ —pV) given by 



(26) 



eM(Ej(k)~p*)/T}-l 



k 2 + rrij and rrij is the j-th meson mass 
Moreover, the boson integrals are sub- 



9j( k ) 

where Ej(k) = 

(see Table HI]) 
jected to the constraint \p*\ < rnj, otherwise Bose con- 
densation becomes possible (as we will see in the next 
Section, this condition is never achieved in the range of 
density and temperature investigated in this paper). The 
values of the effective meson chemical potentials p* are 
obtained from the "bare" ones pj, given in Eq. dT^l) . and 
subsequently expressed in terms of the corresponding ef- 
fective baryon chemical potentials, respecting the strong 
interaction [72]. For example, we have from Eq. (fT5)) that 
Pn+ — p p + — pc = Pp — Pn and the corresponding effec- 



'-P P - PA - (1 - x u]A )g ujN uj - -g pN p , (28) 



Pko(k»°) =Mn _ A*A(£°) 



= Mn - MA - (1 - Xu)A)gu>NU> + -g P NP , (29) 



while the others strangeless neutral mesons have a van- 
ishing chemical potential. Thus, the effective meson 
chemical potentials are coupled with the meson fields re- 
lated to the interaction between baryons. As we will sec 
in the next Section, this assumption represents a crucial 
feature in the EOS at finite density and temperature and 
can be seen somehow in analogy with the hadron res- 
onance gas within the excluded-volume approximation. 
There the hadronic system is still regarded as an ideal 
gas but in the volume reduced by the volume occupied 
by constituents (usually assumed as a phenomenological 
model parameter) , here we have a (quasi free) mesons gas 
but with an effective chemical potential which contains 
the self-consistent interaction of the mesons fields. 

Finally, the total pressure and energy density will be 



P = P B + Pm 
e = e B + e M ■ 



(30) 
(31) 



At a given temperature T, all the above equations must 
be evaluated self-consistently by requiring the baryon, 
electric charge fraction and strangeness numbers conser- 
vation _68|. Therefore, at a given baryon density ps, a 
given Z/A net electric charge fraction (pc = Z/Aps) 
and a zero net strangeness of the system (p$ = 0), the 
chemical potentials pb, pc and ps are univocally deter- 
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TABLE III: Vacuum masses (given in MeV) of the considered hadronic particles. 



M N 


M A 


Me 


M A 


M B 


Mr 


M K 


M v 


M , 
n 




M p 


M u 


M 


939 


1116 


1189 


1232 


1315 


140 


494 


547 


958 


892 


771 


782 


1020 



mined by the following equations 



P [MeV/fm j ] 





= ^2 b i Pi( T ,PB,PC,Ps) , 






GM3 


PB 


(32) 


200 


TM1 




i 




100 


NLp £ 


PC 


= Y2 Cl ft ( T ' PB,PC,PS) , 

i 


(33) 


50 




PS 


= X! St Pt( T 'PB,PC,Ps) , 

i 


(34) 


20 
10 





where the sums run over all considered particles. 



III. RESULTS AND DISCUSSION 

A. Equation of state of pnA symmetric nuclear 
matter 

Let us start our numerical investigation by considering 
the symmetric hadronic EOS using the model discussed 
here. In order to better focalize the role of A-isobar 
degrees of freedom, we will limit this first subsection to 
consider only protons, neutrons and A particles. 

In Fig. [TJ we report the pressure as a function of the 
baryon density (in units of the nuclear saturation density 
po) in the limit of zero temperature. Among the several 
parameters sets, we choose the three sets: GM3, NLp<5, 
TM1 (see Table 1 for details) that meet in a satisfactory 
way the region, reported as shaded area, of pressures con- 
sistent with the experimental flow data of heavy ion colli- 
sions at intermediate energy, analyzed by using the Boltz- 
mann equation model [261 ] . Furthermore, these parame- 
ters sets are largely used to describe the hadronic EOS 
on high density /^-equilibrium compact stars [TH, H3, [HJ . 

At the scope of giving a roughly indication of the pres- 
ence of the A-isobars degrees of freedom from the point 
of view of the stiffness of the EOS, in Fig. [TJ we show 
the behavior of the pressure corresponding only to nu- 
cleons (monotonic curves) and nucleons and A (non- 
monotonic curves) with the scalar r s = g a &/g<jN and 
vector r v = g^A/g^N meson- A couplings ratios. Let us 
note in this context that, when a metastable condition 
for A-isobars is not realized (see discussion below), de- 
cays rates are not taken into account. Moreover, a further 
softening of the EOS can occur considering the degrees of 
freedom of the other hadronic particles. At zero temper- 
ature and symmetric nuclear matter these effects occur 
at very high density and, as anticipated, we will consider 
them separately in the next subsections. 

To better understand the dependence of the EOS on 
the meson-A coupling constants for the different param- 
eters sets, we start by reporting in Fig. [5] the energy per 



r s =\.2, r v =\ 




Pb/Po 



FIG. 1: Pressure as a function of the baryon density (in units 
of the nuclear saturation density po) for the symmetric nuclear 
matter at zero temperature for three different EOSs and A 
coupling ratios (absence of A contribution in the monotonic 
curves). The shaded region corresponds to the limits obtained 
from the analysis of Ref. [2|| . 



baryon versus baryon density at zero temperature and 
GM3 parameters set. The curves (labelled with a, b, c, 
d, e, f) represent different values of the scalar r s and 
vector r v meson-A couplings ratios. In setting these cou- 
pling constants we have required, as in Ref.[46j, that i) 
the second minimum of the energy per baryon lies above 
the saturation energy of normal nuclear matter, i.e., in 
the mixed A-nucleon phase only a metastable state can 
occur; ii) there are no A-isobars present at the satura- 
tion density; iii) the scalar field is more (equal) attractive 
and the vector potential is less (equal) repulsive for As 
than for nucleons, in according with QCD finite density 
calculations [Hj]. Of course, the choice of couplings that 
satisfies the above conditions is not unique but exists a 
finite range of possible values (represented as a triangle 
region in the plane r s -r v ) which depends on the partic- 



ular EOS under consideration [46j. Without loss of gen- 
erality, we can limit our investigation to move only in a 
side of such a triangle region by fixing r v = 1 and vary- 
ing r s from unity to a maximum value compatible with 
the conditions mentioned above. Similar conclusions are 
obtained with any other compatible choice of the two 
coupling ratios. 

In Fig. [3J we compare the energy per baryon for dif- 
ferent parameters sets but at a fixed value of the scalar 
and vector A coupling constants. At variance of the pa- 
rameters sets we have a very different behavior, on the 
other hand, we obtain comparable features for the three 
considered parameters sets by means of a rescaling of 
the A couplings. To better clarify this aspect, in Table 
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Pb/Po 



FIG. 2: The energy per baryon versus baryon density at zero 
temperature and GM3 parameters set with: a) without A, 
b) non-interacting A (r s — r v = 0), c:(r 3 = 1.3, r v = 1), 
d:{r a = 1.41, r v = 1), e:(r s = 1.45, r v = 1), f:(r s = 1.5, 
r„ = 1). 

e/ps-M N [MeV] 




FIG. 3: (Color online) The same as Fig. [2] but for different 
parameters sets and fixed r s — 1.3 and r v = 1. 



IIVI we report, for the three parameters sets and fixing 
r v = 1, the values r' 1 corresponding to the appearance 
of the second minimum on the energy per baryon and 
r™ corresponding to the maximum values of r s com- 
patible with the constraint that the second minimum of 
the energy per baryon lies above the saturation energy of 
normal nuclear matter. 

In order to get a deeper insight on the dependence 
of the A-isobars from the coupling constants and from 
the temperature, in Fig. 01 we show the behavior of 
the effective mass for different temperatures and A 



TABLE IV: Values of the r s corresponding to the appearance 
of the second minimum on the energy per baryon and the 
maximum values r™ ax obtained by requiring that in the mixed 
A-nucleon phase only a metastable state can occur. 





GM3 


NLp<5 


TM1 


' S 


1.41 


1.32 


1.27 


.max 
' s 


1.50 


1.41 


1.33 



M* A /M A 




FIG. 4: The A effective mass ratios versus the baryon density 
with r s = r v = 1 (upper curves) and with r s — 1.32 and r v = 
1 (lower curves) for different temperatures. The parameters 
set used in the calculation is NLp<5. 



coupling constants (r s = r v = 1 for upper curves and 
r s = 1.32 and r v = 1 for lower curves) in the NLp<5 
parameters sets. Let us note that the coupling ratio of the 
lower curves corresponds to the appearance of the second 
minimum of the energy per baryon at zero temperature 
(metastable state). As we can see, the different behavior 
of the effective A mass from the coupling constants is 
strongly linked to the different behavior of the energy 
per baryon versus baryon density and temperature. 

In Fig. [5J we show the energy per baryon as a func- 
tion of the baryon density at T = 40 4- 100 MeV for two 
different A scalar coupling ratios r s and for the TM1 
parameters set. As stated before, at zero temperature 
and fixed r v = 1, there is no second minimum if r s = 1 
while it takes place if r s = 1.3 (see Fig. [3] and Table HV)l . 
As the temperature increases, the behavior is remarkable 
different in the case of r s = 1.3 where the first minimum 
disappears and only the second minimum remains. In the 
above two figures, the results are reported only for one 
parameters set, however similar behaviors are obtained 
with the other sets, if we use comparable meson-A cou- 
pling ratios. To better clarify this aspect, in Fig. [51 we 
report the variation of the baryon density, with respect to 
temperature, corresponding to the position of the second 
minimum of the energy per baryon for the three different 
parameters sets. In the comparison, we fix r s to the av- 
erage value r m between r^ 1 and r™ ax listed in Table IIVI 
for each parameters sets (r s = 1.46 for GM3, r s — 1.37 
for NLp(5, r s = 1.30 for TM1). For all three sets, the 
position of the second minimum appears approximately 
at a constant value of baryon density ps ~ 3 -j- 4 po until 
T « 80 -j- 100 MeV. At higher temperatures, as observed 
in Fig. [5l the first minimum disappears and the second 
minimum moves rapidly at lower baryon densities. At 
fixed temperature, the different positions of the second 
minimum are a direct consequence of the different sat- 
urated nucleon effective mass Mjy- in the considered pa- 
rameters sets. In agreement with the results of Ref . [441] . 
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e/p B -M N [MeV] 




FIG. 5: The energy per baryon versus baryon density for 
symmetric nuclear matter at different values of temperature 
and TM1 parameters set. The solid lines correspond to the A 
coupling ratios r 3 — r v = 1 and the dashes lines correspond 
to r s — 1.3 and r v = 1. 



Pb/Po 




FIG. 6: (Color online) Variation of the baryon density, with 
respect to temperature, corresponding to the position of the 
second minimum of the energy per baryon for different pa- 
rameters sets (see text for details). 



we have, in fact, that a smaller value of favors the ap- 
pearance of a second minimum of the energy per baryon 
at lower baryon density. Therefore, the position of the 
second minimum results to be at a sensibly lower baryon 
density for the TM1 set with respect to the NL/?<5 and 
GM3 ones. 

In Fig. [7J we report the nucleon (solid lines) and the 
A-isobar (dashed lines) density, normalized to the baryon 
density, versus the baryon density for different values of 
temperature (in unit of MeV). We observe that, although 
the A-isobar density seems to be negligible at low tem- 
peratures up to very high densities, it becomes relevant 
by increasing the temperature even for r s = r v = 1 (up- 
per panel). Moreover, in the case of r s = 1.32 and r v = 1 
(lower panel), the A particle density becomes compara- 
ble to the nucleon density in the range of T w 80 -j- 120 
MeV and ps ~ 1 -r- 2. 5 po; values that can be reached in 
high energy heavy ion collisions. This behavior has been 
obtained for the NLpS parameters set but is common for 
all three considered sets even if, at fixed temperature and 



pIpb 




FIG. 7: The relative nucleon (solid lines) and A (dashed lines) 
density versus the baryon density with r s — r v = 1 (upper 
panel) and with r s = 1.32 and r v = 1 (lower panel) for differ- 
ent values of temperature (in units of MeV). The parameters 
set used in the calculation is NLp5. 



baryon density, different values of particle densities are 
obtained for different EOSs. To better focus this matter 
of fact, in Fig. [8j we report the variation of the baryon 
density, as a function of temperature, for which A-isobar 
density is equal to nucleon density (pa = Pn = ps/2), 
for the three different parameters sets. Also in this case, 
in the comparison, we use comparable values of r s for 
different EOSs. As in the previous figure, for the NLpS 
parameters set we fix r s = 1.32, corresponding to the 
value r* 1 of Table HVl Analogously, we use r s = 1.41 for 
GM3 and r 8 = 1.27 for TM1. As already observed in Fig. 
[51 for the TM1 parameters set, which has a lower value of 
M^, the formation of A-isobars occurs at lower baryon 
densities with respect to NLpS and GM3 (with larger val- 
ues of Mjy-). Small variations between NLpS and GM3 
correspond mainly to the almost equal saturated nucleon 
effective mass (slightly greater for the GM3 parameters 
set). 

Finally, in agreement with previous investigations 
[43l [il], let us observe that, in Fig. [71 the pa/pb and 
Pn / Pb ratios become constant at sufficiently high baryon 
density, regardless of the temperature. Increasing the 
r a ratio (and fixing r v ) such constant asymptotic values 
are reached at lower baryon density. Moreover, we have 
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FIG. 8: (Color online) Variation of baryon density, as a func- 
tion of temperature, for which A-isobar density results to be 
equal to nucleon density, for different parameters sets. The 
scalar coupling ratios for each EOS are the values r s — rl 1 
reported in Table HVl 

verified that this behavior is still realized in asymmetric 
hadronic matter and even in the presence of hyperons 
and mesons degrees of freedom. This feature could be 
an interesting matter of investigation in the future high 
energy compressed nuclear matter experiments. 



/UsIMeV] 




B. Equation of state with strange particles 

Let us now investigate the EOS with the inclusion of 
hyperons, non-strange and strange mesons particles at 
fixed values of Z/A and zero net strangeness, as described 
in Sec. II. 

In Fig. |9l we report the isotherms of the strange chem- 
ical potential /is (upper panel) and the electric charge 
chemical potential pc (lower panel) for different values 
of temperature and Z/A = 0.4. To point out the role 
of the As degrees of freedom, we have considered three 
different cases: i) the solid lines do not contain A con- 
tribution, ii) in the long dashed lines the A couplings 
are r s = r v = 1, iii) r s = 1.3 and r v — 1 for the short 
dashed lines. As expected, for a multi-composed strange 
hadronic matter, ps is positive and increases with T and 
fiB- At low T we observe very small variations in the 
strangeness chemical potential with different A coupling 
constants. Very significant effects instead are present in 
the behavior of pc, where, in the presence of the A de- 
grees of freedom, we have a sensible reduction of its abso- 
lute value and remains almost constant at high /ib- This 
matter of fact suppresses the possibility of pion conden- 
sation also at very high baryon density (see below for a 
further discussion on Bose condensation). 

Analogously, in Fig.slTOlandfTTl we report, for different 
temperatures and Z/A = 0.4, the baryon density and the 
pressure as a function of the baryon chemical potential. 
We can see how the presence of the A degrees of freedom 
becomes, already for T « 80 MeV, very remarkable for 
baryon densities greater than about po/2. The parame- 



FIG. 9: Variations of the strangeness chemical potential /is 
(upper panel) and electric charge chemical potential pc (lower 
panel) with respect to the baryon chemical /is at different val- 
ues of temperature (in units of MeV) and different A coupling 
constants. The parameters set is GM3 and r v = 1. 



ters set used in the above calculations is GM3, however 
very similar behaviors can be obtained with the other 
sets. As already remarked, the most relevant difference 
is that the presence of A particles occurs at lower baryon 
chemical potentials (baryon densities) for the TM1 EOS. 

To better understand the relevance of the A-isobars, 
together with the effects of the electric charge fraction 
and of the effective meson chemical potentials, in Fig. 
[T21 we report the relative difference of the pressure AP/P 
(without and with As contribution) as a function of the 
baryon chemical potential. The used parameters set is 
NLp<5 but we have common behaviors for all three sets. In 
the panel (a) and (b) we show the sensibility of the EOS 
with respect to a variation of Z/A. In fact in the panel 
(a) AP(Z/A,n*) = P(Z/A = 0.5) - P(Z/A = 0.4) and 
the symbol /i* means that we have taken account of effec- 
tive meson chemical potentials, as described in Sec. II. In 
the panel (b) is reported the same relative difference but 
with bare meson chemical potentials /i (in other words, 
all mesons are considered as a free gas of non-interacting 
particles). In the figure, the circles, the squares and the 
triangles represent the values of the baryon chemical po- 
tential corresponding to ps = Poj 2po and 3po, respec- 
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FIG. 10: Baryon density (in units of the nuclear saturation 
density po) as a function of the baryon chemical potential 
/ib at different values of temperature (in units of MeV). The 
parameters set is GM3 and r v — 1. 



P [MeV/fm j ] 




FIG. 11: Pressure as a function of the baryon chemical poten- 
tial /is at different values of temperature (in units of MeV). 
The parameter set is GM3 and r v = 1. 



tively. As expected, the EOS is more sensible to a varia- 
tion of Z/A at low temperature and this effect decreases 
by increasing the temperature. However, this behavior is 
strongly related to the A-isobars degrees of freedom. The 
presence of A-isobars greatly reduces the dependence on 
Z/A in the range of baryon density and temperature rel- 
evant in our investigation. This matter of fact is also 
realized by considering bare meson chemical potentials 
(panel (b)), even if this effect is less marked. As a con- 
sequence, we expect that A-isobars degrees of freedom 
affect significantly the value of the symmetric energy at 
finite density and temperature. 

In the panel (c) and (d) of Fig. Q21 we wish to empha- 
size the importance of effective meson chemical potentials 
by considering the relative difference between the pres- 
sure P(p), calculated with bare meson chemical poten- 
tials, and the pressure P(p*), including effective meson 
chemical potentials at fixed ratio Z/A = 0.4 (panel (c)) 
and Z/A — 0.5 (panel (d)). The presence of effective me- 
son chemical potentials reflects the behavior of the self- 
consistent values of the meson fields and, in particular, we 



have that its relevance: i) depends on the baryon density 
(or /is), ii) increases with the temperature, iii) depends 
on the isospin asymmetry (more relevant for asymmet- 
ric hadronic matter), iv) decreases with the presence of 
A and the relative difference has a maximum in the re- 
gion of pb ~0.5-^3po- We will see in next subsection 
as these features will be very important in the behavior 
of the considered particle-antiparticle ratios and in the 
strangeness production. 

Always concerning the role of the effective meson 
chemical potential, let us further observe that its abso- 
lute value is significantly lower than its bare value and, 
therefore, the window of /is and T values, in which Bosc 
condensation can occur, appears to be greatly reduced. 
Considering, for example, K + meson, its bare chemical 
potential (see Eq.([I5|) is p K + = fj,g + p c = p p — p\, 
which is dominated from the behavior of ps (that in- 
creases with [Lb and T, see Fig. However, taking into 
account Eq.(|28p. [i* K+ is significantly reduced compared 
to fi K + by the presence of the uj and p mesons fields. 
At fixed ratio Z/A = 0.4 and zero net strangeness, by in- 
creasing pb and T, the term containing the p meson field 
is always negative but is much smaller than the (positive) 
u) meson field one and, thus, kaon condensation can occur 
only at very high baryon densities Q . 



C. Particle ratios 

We start this subsection by considering, in Fig. [T3l 
the ratio of the net densities A ++ /p as a function of 
the baryon chemical potential for different values of tem- 
perature and different parameters sets. As before, we 
fix Z/A = 0.4. Let us observe that the ratio increases 
with the temperature but is almost constant with ps for 
r s = 1. On the other hand, it increases rapidly with ps if 
we increase the A coupling r s . In agreement with the pre- 
vious results of Fig. s [6] and HI the formation of A particles 
appears to be strongly favored for the TM1 parameters 
set with respect the other two, especially for r s = 1.3. 
Let us remember that, in this last case, A particles are in 
different regime depending on the considered parameters 
sets (see Table HV|). The above behavior should be espe- 
cially evident at low transverse momentum pion spectra 
in heavy ion collisions at intermediate/high baryon den- 
sities and temperatures. Furthermore, remembering that 
we are not considering decays and rescattering effects, we 
observe that the order of magnitude of the ratio at high 
temperature and low density seems to be compatible with 
the SPS/RHIC experimental results (with measured ra- 
tios approximately equal to 0.2 -j- 0.5) [HI |4^. 

In Fig. Q3J we show the variation of K + /tt + and 
K~/tt~ ratios with respect to temperature at various 
baryon chemical potentials, considering different parame- 
ters sets at fixed r s = r v = 1 coupling ratios. Appreciable 
variations between the EOSs are observed only at higher 
baryon chemical potentials (pb — 600 MeV). By increas- 
ing ps we have that the difference between K + /ir + and 
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FIG. 12: Relative difference of the pressure as a function of the baryon chemical potential ps at different values of temperature 
with the exclusion (no A) and the inclusion (r s = 1.3 and r v = 1) of the A-isobars degrees of freedom. Panel (a): AP(Z/A, fj,*) = 
P(Z/A = 0.5) — P(Z/A = 0.4) and the symbol p* means that we have considered an effective chemical potential for all hadrons 
(see text for details); panel (b): the same of (a) but mesons have a bare chemical potential (free boson gas); panel (c): 
AP(Z/A = 0.4, p, /i*) = P(p) — P(p*) at fixed Z/A=0.4, where the pressure P(p) is calculated by considering a bare chemical 
potential for mesons and the pressure P(p*) is calculated by using an effective chemical potential for all hadrons; panel (d): 
same of (c) but at fixed Z/A=0.5. The circles, the squares and the triangles represent the values of the baryon chemical 
potential corresponding to pB — po, 2po and 3po, respectively. 



K~ /ir~ ratios increases with the temperature but such 
a difference becomes much smaller at low hb- This be- 
havior is in agreement with recent relativistic heavy ion 
collisions data [Til ] . 

In Fig. [T5l we show the ratio K + /K~ as a function of 
temperature at different ps and A coupling ratios (solid 
lines: absence of As, long dashed lines: r s — r v = 1, 
short dashed lines: r s — 1.3, r v = 1); Z/A = 0.4. As 
expected, we have a value of the ratio nearly equal to 
one at low baryon chemical potentials (pis < 300 MeV), 
while, for higher ps, the ratio has a peak corresponding 
to baryon density ps ~ 0.1 -f- 0.2 p for ps ~ 500 -=- 
600 MeV curves. This non-monotonic behavior is much 
more evident taking into account the A-isobars degrees 
of freedom. 

In order to investigate how the previous results depend 
on the choice of the EOS, in Fig. (THl we report the 
variation of the K + /K~ ratio with respect to baryon 
chemical potential, at fixed temperature T = 100 MeV 



and for the parameters sets GM3 and TM1 (the results 
relative to the NLpS set are not reported because very 
close to the GM3 ones). As in Fig. [51 we fix the higher 
value of r s to the average value r m between r* 1 and r™ ax 
listed in Table[TV](r m = 1.46 for GM3 and r m = 1.30 for 
TM1 parameters set). The two EOSs have similar shapes 
even if the ratio results to be suppressed at higher /is in 
the TM1 model, also in absence of A degrees of freedom. 
Moreover, a very peculiar behavior appears at r s = r m 
where the ratio has a peak around ps ~ 640 -j- 650 MeV 
(corresponding to pb ~ 0.7po and ps ~ Po for TM1 and 
GM3 sets, respectively). This feature is mainly due to 
the fact that, at fixed Z/A and zero net strangeness, the 
density of K + mesons are reduced by the presence of A 
particles. 

Related to the above results, it should be interesting to 
investigate from the experimental point of view the be- 
havior of the ratio K + / K~ in a kinematical region cor- 
responding to intermediate/high temperatures and high 
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FIG. 13: (Color online) Ratios of net densities A ++ /p as a 
function of the baryon chemical potential for different temper- 
atures, r a and parameters sets (a: TM1, b: NLpS, c: GM3). 



FIG. 16: (Color online) Variation of the K + /K~ ratio with 
respect to baryon chemical potential at fixed temperature 
T — 100 MeV. The value r m corresponds to the average value 
between r\ l and r" lax listed in Table [IV] for the two different 
parameters sets. 



K/tt 




FIG. 14: (Color online) Variation of the K + /ir + and K /ty 
ratios with respect to temperature at different values of 
baryon chemical potential and different parameters sets (a: 
TM1, b: NLpS, c: GM3). The A coupling ratios are fixed to 
r a = r v = 1. 
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FIG. 15: Variation of the K + / K~ ratio with respect to tem- 
perature at different values of baryon chemical potential and 
for different A coupling ratios (TM1 parameters set). 



values oi hb- 

In order to gain a deeper insight about the role of the 
A-isobars and the net electric charge fraction, we report 
in Fig.s[l2]and[T8]the variation of K+/tt+ and K + /K~ 
ratios with respect to baryon density at T = 80 and 
120 MeV. For each temperature, we have set: Z/A=0A 
for solid curves and Z/A=0.5 for long dashed curves, 
both in absence of A particles; Z/A=0A and A cou- 
plings r s = 1.2 for short dashed curves; Z/A=0.5 and 
r s = 1.2 for dotted curves (r v = 1). The parameters set 
is TM1 but similar behaviors are observed for the other 
two parameters sets. 

Concerning the dependence on Z/A, we have to com- 
pare solid with long dashed curves (in absence of A par- 
ticles) and short dashed with dotted ones (with the in- 
clusion of As). We can see that the differences between 
K + /ir + ratios are comparable for the two temperatures, 
while we have a decreasing difference between K + /K~ 
ratios by increasing the temperature (it occurs an en- 
hancement of K + / K~ ratios at Z/A = 0.5 with respect 
to the value Z/A = 0.4; the other way round takes place 
for K + /tt + ratios). As already observed in Fig. [12] the 
presence of the A-isobars strongly suppresses the depen- 
dence on Z/A for the considered particle ratios. On the 
other hand, concerning the strangeness production in the 
presence of the A-isobars degrees of freedom, we can com- 
pare solid with short dashed curves (at fixed Z/A = 0.4) 
and long dashed curves with the dotted ones (at fixed 
Z/A — 0.5). In agreement with the results of Fig.s [T4l 
[TBI and [TBI we observe that, at fixed T and ps, the pres- 
ence of the A-isobars sensibly decreases both K + /ir + and 
K + /K~ ratios. Moreover, to outline the importance of 
the effective meson chemical potential fi*, in Fig.sfTTland 
[TBI we have inserted the dash-dotted curves correspond- 
ing to ratios with Z/A=0A and A coupling r s = 1.2, 
but with bare meson chemical potentials /i. As already 
outlined in Fig. 1121 comparing these last curves with the 
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short dashed ones, it is possible to observe the relevance 
of the effective meson chemical potentials at finite den- 
sity and temperature and how they thus avoid unphysical 
too high ratios (7l| . 

Finally, it is interesting to extend the study of the 
EOS also at high temperatures and low baryon chemi- 
cal potentials regime. At this scope, in Fig. [19J we re- 
port the results of various particle-antiparticle ratios and 
K + /ir + ratio as a function of the p/p ratio for different 
values of temperature. The A coupling ratios are fixed 
to r s = r v = 1 and Zj A = 0.4. The ratios are reported 
for the GM3 parameters set, however, we have verified 
that very close results are obtained for the other two pa- 
rameters sets. Also in this case we can observe good 
agreements with the results obtained in the framework 
of statistical-thermal models 18[ and with experimental 
SPS and RHIC data [hI 
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FIG. 17: Variation of K + /n + ratios with respect to baryon 
density at T = 80 MeV (lower curves) and T = 120 MeV 
(upper curves). For the dash-dotted curves the symbol fi 
indicates that all mesons have bare chemical potentials (see 
text for details). 
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FIG. 18: The same of Fig[l7]but for K+/R- ratios 
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FIG. 19: Particle-antiparticle and K + /tv + ratios as a function 
of the p/p ratio for different temperatures. The A coupling 
ratios are fixed to r s = r v — 1. The ratios of E + /H~ at 
T = 80, 120 MeV are not reported because very strictly to 
the A/A ones. 



IV. CONCLUSIONS 

The main goal of this paper is to show systematically 
how the presence of the A-isobars degrees of freedom 
affects the hadronic EOS by requiring, in the range of 
finite density and temperature, the global conservation 
of baryon number, electric charge fraction and zero net 
strangeness. In this study we have considered three dif- 
ferent parameters sets (GM3, NLp<5, TM1) that are com- 
patible with recent analysis at intermediate energy heavy 
ion collisions and extensively adopted in several applica- 
tions related to high density /3-equilibrium compact stars. 
We have studied a RMF model with the inclusion of the 
full octet of baryons and A-isobars, self-interacting by 
means of a, p, u>, 5 mesons fields. To take into account 
the lightest pseudoscalar and vector mesons contribution, 
especially in regime of low (but finite) baryon density 
and high temperature, we have incorporated the mesons 
as an ideal Bose gas but with effective chemical poten- 
tials. We have shown that in the EOS and, as a conse- 
quence, in the analyzed particle ratios, this assumption 
appears to be very relevant in the range of density and 
temperature considered in this paper. The role of the 
effective meson chemical potential has a phenomenolog- 
ical counterpart in the excluded volume approximation 
for the hadron resonance gas where all effective particle 
chemical potentials are shifted, respect to the real ones, 
proportionally to the particle excluded volume (usually 
assumed as a parameter). Here, from a more microscopic 
point of view, the effective meson chemical potentials are 
shifted proportionally to the mesons fields, related to the 
self-consistent interaction between baryons. 

The relevance of the A-isobars in the EOS has been 
investigated for different parameters sets and coupling 
ratios at zero and finite temperature, in absence and 
in presence of hyperons and mesons, for symmetric and 
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asymmetric nuclear matter. In all considered cases we 
have shown that the A-isobars degrees of freedom play 
a crucial role. Especially in the range of finite density 
and temperature considered in the last two subsections, 
we have seen that: i) the relevance of A-isobars strongly 
increases with a density and temperature increase; ii) at 
fixed Z/A, the presence of A-isobars in the EOS affects 
significantly the strangeness production; iii) A-isobars 
degrees of freedom remarkably decrease the dependence 
on the isospin of the EOS and, as a consequence, of the 
considered particle ratios. This last property appears to 
be very relevant also in connection to the supposed de- 
pendence on Z/A of the critical transition density from 
hadronic matter to a mixed phase of quarks and hadrons 
at high baryon and isospin density. 

All quoted features are realized even if the A-coupling 
ratios r s and r v do not necessarily correspond to the for- 
mation of a A-isobar metastable state, however, much 
stronger effects are present if we consider r s > r), 1 indi- 
cated in Table IPV1 Furthermore, apart from an appropri- 
ate variation of coupling constants, the obtained results 
are comparable for all three considered parameters sets. 



The most relevant difference is that the formation of A 
particles is favored at lower baryon density for the TM1 
EOS with respect to the other two parameters sets. This 
matter of fact is strongly correlated to the significantly 
lower value of the saturated effective nuclear mass 
in the TM1 parameters set. 

Finally, we have investigated several particle- 
antiparticle ratios and strangeness production as a func- 
tion of the temperature, baryon density and antiproton 
to proton ratio. Although the studied EOS is princi- 
pally devoted to a regime of finite and intermediate val- 
ues of baryon density and temperature, a satisfactory 
agreement has been found with recent relativistic heavy 
ion collisions data and with statistical-thermal models 
results. 
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